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ABSTRACT 


The  Special  Sensor  Microwave/imager  (SSM/l)  will  be  launched 
aboard  a  DMSP  satellite  in  FY-85.  Seven  channels  of  passive 
microwave  data  will  be  provided  by  the  sensor  at  frequencies  of 
19.35  GHz,  22.235  GHz,  37.0  GHz,  and  85.5  GHz.  Dual 
polarizations  will  be  measured  in  all  but  the  22.235  GHz  channel 
where  only  the  vertically  polarized  radiation  will  be  sensed. 
Pixel  size  is  dependent  on  frequency  and  ranges  from  approximate¬ 
ly  55  km  square  at  19.35  GHz  to  less  than  15  km  square  at  85.5 
GHz.  Spatial  resolution  of  the  high  frequency  85.5  GHz  channel 
is  substantially  better  than  the  highest  resolution  obtained  from 
ESMR  and  SMMR  sensors  (approximately  30  km).  Data  received  from 
this  channel  will  improve  the  accuracy  with  which  ice  edges  and 
bodies  of  open  water  within  the  pack  can  be  located.  The  cell 
size  for  which  total  ice  concentration  can  be  estimated  also  will 
be  reduced.  However,  estimates  of  multi-year  ice  concentration, 
concentration  of  thin  ice  types,  and  the  character  of  ice  in  the 
marginal  ice  zone  will  not  be  improved  over  SMMR  estimates  by 
virtue  of  these  higher  frequency  data.  Enhanced  estimates  of 
these  parameters  must  come  from  development  of  improved  ice 
classification  algorithms.  Improved  ground  truth  information 
from  which  satellite  data  can  be  interpreted  and  verified  is 
critical  if  classification  algorithms  are  to  be  improved. 
Specifically,  data  that  document  radiometric  signatures  of  wet 
and  snow-covered  ice  surfaces  that  occur  during  summer  and  autumn 
months  are  needed  to  improve  estimates  of  multi-year 
concentrations  during  melt  and  freeze-up.  Brightness  temperature 
data  that  describe  the  transition  from  open  water  to  thick 
first-year  ice  are  required  for  improved  discrimination  of  thin 
ice  types.  Ground  truth  data  that  correlate  marginal  ice  zone 
conditions  with  low-re so  1 ut i orj  imagery  from  satellite-borne 
microwave  sensors  are  needed  to  improve  estimates  of  conditions 
within  the  marginal  ice  zone. 
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INTRODUCTION 


Data  obtained  from  satellite-borne  imaging  microwave  radiome¬ 
ters  are  relevant  to  academic  and  applied  problems  in  polar  re¬ 
gions.  As  all-weather  day-night  sensors,  they  provide  coverage 
of  Arctic  and  Antarctic  regions  regardless  of  surface  conditions. 
Information  concerning  sea  ice  concentration  and  age,  daily  con¬ 
figuration  of  the  ice  edge,  location  of  polynyas,  movement  of  ice 
masses  and  weather  systems,  and  estimates  of  heat  exchange 
between  ocean  and  atmosphere  in  polar  regions  is  derived  from 
these  data. 

Launch  of  the  Special  Sensor  Microwave/imager  (SSM/l)  in 
1985  aboard  the  DMSP  satellite  will  provide  a  new  source  of 
all-weather,  all-light  images  of  the  polar  regions.  Projected 
launches  of  consecutive  SSM/l'  s  on  DMSP  platforms  and  the  Navy 
Remote  Ocean  Sensing  System  (NROSS)  platform  at  the  end  of  the 
decade  present  the  possibility  of  nearly  continuous  polar  cover¬ 
age  for  the  next  ten  years.  Potential  benefits  accrued  from 
these  sensors  will  be  reaped  in  the  operational,  tactical, 
academic,  and  research  communities.  Development  of  new  analysis 
procedures  that  result  could  lead  to  significant  improvements  in 
estimating  ice  conditions  in  the  marginal  ice  zone,  more  accurate 
predictions  of  multi-year  ice  concentrations  during  summer  and 
autumn  months,  better  assessment  of  heat  balance  parameters  in 
ice  covered  regions,  and  a  significantly  improved  understanding 
of  basin-scale  dynamics  of  the  Arctic  pack. 


HISTORICAL  DEVELOPMENT :  ESMR-SMMR-SSM/l 

The  Electrically  Scanning  Microwave  Radiometer  (ESMR-5)  was 
launched  aboard  Nimbus-5  in  December  1972  to  provide  maps  of 
precipitating  clouds  over  ocean  areas  and  to  chart  the  distribu¬ 
tion  of  sea  ice  in  polar  regions  (Sabatini,  1972).  ESMR  utilized 
a  single  frequency  (19.35  GHz)  passive  microwave  radiometer  with 
a  beam  width  of  1.4°  x  1.4°  at  nadir.  The  instrument  employed  a 
phased  array  antenna  to  scan  50°  either  side  of  the  sub-satellite 
point.  A  2500  km  wide  swath  was  included  in  each  scan. 
Microwave  images  returned  by  the  sensor  are  composed  of  bright¬ 
ness  temperatures  measured  at  each  of  78  positions  across  the 
scan.  Spatial  resolution  of  these  images  is  25  km  x  25  km  at 
nadir,  and  degrades  to  160  km  crosstrack  x  45  km  downtrack  at  the 
ends  of  the  scan.  Polar  coverage  was  virtually  continuous  due  to 
the  wide  swath  width;  overlap  between  successive  passes  covered 
approximately  50%  of  latitudes  above  70°.  ESMR-5  functioned  well 
until  1976  when  certain  beam  positions  in  the  phased  array 
antenna  were  lost.  As  late  as  1980  ESMR-5  was  the  primary  source 
of  data  used  to  compile  Navy  operational  ice  analyses  (Wilheit, 
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1980).  An  expanded  discussion  of  the  ESMR  instrument  is  given  by 
Wilheit  (1972) . 

During  its  four  year  effective  life,  ESMR-5  returned  data 
that  established  the  value  of  satellite-borne  microwave  imagers 
to  polar  studies.  Numerous  workers  have  shown  that  diverse 
information  concerning  total  sea  ice  concentration  (Gloersen,  et 
al . ,  1974,  1978;  Campbell,  et  al * ,  1980;  Comiso  and  Zwally,  1982; 
Carsey,  1982),  multi-year  ice  concentration  (Carsey,  1982),  con¬ 
figuration  of  the  ice  edge  (Gloersen,  et  al . ,  1974;  Carsey, 
1982),  location  of  polynyas  (Zwally  and  Gloersen,  1977;  Rayner 
and  Howarth,  1979;  Carsey,  1980),  and  changes  in  the  structure 
and  character  of  the  Arctic  ice  pack  with  time  (Campbell,  et  al . , 
1980)  is  contained  in  ESMR-5  imagery.  ESMR-5 1 s  primary  weakness 
was  its  single-channel  capability.  Real-time  discrimination  of 
first-year  ice  from  multi-year  ice  requires  that  radiation  at 
several  frequencies  be  compared.  Furthermore,  lack  of  accurate 
information  regarding  atmospheric  effects  that  can  be  obtained 
from  multi-channel  radiometers  limits  the  precision  with  which 
concentration  estimates  can  be  determined  from  the  single-channel 
ESMR  sensor. 

ESMR-6,  a  successor  to  ESMR-5,  was  launched  aboard  Nimbus-6 
in  June  197  5.  The  ESMR-6  instrument  was  designed  as  a 
multi-channel  sensor  and  measured  brightness  temperatures  at  37 
GHz  in  both  horizontal  and  vertical  polarizations.  Failure  of 
the  horizontal  channel  in  September  1976  reduced  ESMR-6  to  a 
single-channel  sensor.  Although  data  obtained  by  ESMR-6  are 
potentially  useful  for  mapping  the  extent  of  ice  coverage  and 
discriminating  between  first-year  and  multi-year  ice  types, 
ESMR-6  data  largely  have  not  been  used  in  polar  studies. 

The  Scanning  Multichannel  Microwave  Radiometer  (SMMR)  flown 
aboard  the  Nimbus-7  (Nimbus-G)  and  Seasat-A  satellites  in  1978 
was  designed  to  rectify  weaknesses  inherent  in  the  ESMR  instru¬ 
ment.  Technology  associated  with  the  Nimbus  ESMR,  Microwave 
Spectrometer  (NEMS),  and  Scanning  Microwave  Spectrometer  (SCAMS) 
were  incorporated  into  SMMR  design  (Gloersen  and  Barath,  1977  ). 
The  SMMR  instrument  measures  microwave  radiation  through  ten 
channels  at  five  frequencies.  Passive  radiation  is  sensed  in 
both  vertical  and  horizonatal  polarizations  at  6.633,  10.69, 
18.0,  21.0,  and  37.0  GHz.  For  the  37  GHz  channels,  two  radiome¬ 
ters  are  used,  one  for  each  polarization.  At  each  of  the  other 
four  frequencies  a  single  radiometer  alternates  between  polariza¬ 
tions  on  sucessive  scans.  The  angle  of  view  for  the  instrument 
is  42°  which,  in  the  case  of  Seasat-A,  produced  a  swath  width  of 
595  km.  Cell  dimensions  for  a  single  pixel  depend  upon  channel 
frequency  and  vary  from  149  km  x  149  km  at  6.6  Ghz  to  an  optimum 
27  km  x  27  km  at  37  GHz  for  Seasat-A.  The  instrument  aboard 
Nimbus-7  is  forward  viewing  with  the  swath  centered  on  the 
subsatellite  track.  The  Seasat-A  instrument  was  aft  viewing  with 
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the  swath  center  22°  from  nadir  to  the  right  of  the  satellite  to 
permit  overlap  with  footprints  of  the  onboard  scatterometer  and 
synthetic  aperture  radar.  Both  Nimbus-7  and  Seasat-A  were 
launched  in  1978.  Seasat-A  operated  for  95  days  from  July  6  to 
October  10.  Operation  ceased  when  the  onboard  power  system 
failed.  Nimbus-7  has  functioned  from  1978  to  present.  Technical 
discussions  of  SMMR  characteristics  are  presented  by  Gloersen  and 
Barath  (1977),  Dunne  (1978),  Barrick  and  Swift  (1930),  Njoku,  et 
al .  (1980a,  1980b),  and  by  Swanson  and  Riley  (1980). 

Improved  spectral  resolution  offered  by  the  multi-channel 
SMMR  instruments  permitted  development  of  algorithms  that 
significantly  enhance  ice  concentration  prediction  capabilities 
(Gloersen,  e t  a  1 .  ,  1981a,  1981b;  Svendsen,  e t  a  1 . ,  1983  ). 
Real-time  multi-year  ice  concentration,  a  quantity  that  could  not 
be  calculated  daily  from  ESMR  data,  is  calculated  to  an  accuracy 
of  +_10%  from  SMMR  data.  Total  ice  concentration  estimates  are 
accurate  to  +3%  under  ideal  conditions  (Svendsen,  et  al . ,  1983). 
Location  of  the  ice  edge  can  be  estimated  to  within  HyLO  km. 
However,  detail  within  the  marginal  ice  zone  can  not  be  resolved 
due  to  the  large  30  km  square  footprint  of  the  highest  frequency 
channel.  For  similar  reasons,  many  tactically  significant  bodies 
of  open  water  can  not  be  delineated. 

The  Special  Sensor  Microwave/imager  (SSM/l)  will  be  launched 
aboard  the  DMSP  satellite  in  FY-85.  Seven  channels  of  passive 
microwave  data  at  frequencies  of  19.35  GHz,  22.235  GHz,  37.0  GHz, 
and  85.5  GHz  will  be  provided.  Both  vertical  and  horizonal  po¬ 
larizations  will  be  measured  in  all  but  the  22.235  GHz  channel 
where  only  vertically  polarized  data  will  be  sensed.  The 
effective  footprint  size  ranges  from  55  km  for  the  low  frequency 
19.35  GHz  channel  to  15  km  for  the  high  frequency  85.5  GHz 
channel.  Spatial  resolution  provided  by  the  85.5  GHz  channel  is 
improved  significantly  over  SMMR  and  should  provide  refined 
estimates  of  ice  edge  locations.  Polynyas  of  smaller  size  than 
previously  seen  in  ESMR  and  SMMR  imagery  should  be  seen  in  SSM/l 
images.  Accuracy  of  ice  concentration  estimates  is  projected  by 
Hughes,  Inc.  to  be  +12%.  Ice  edges  are  expected  to  be  located  to 
+12.5  km.  In  light  of  accuracies  achieved  with  SMMR  data,  SSM/l 
estimates  should  exceed  these  projections,  especially  after 
algorithms  are  tuned  with  ground  truth  data. 


ICE  CONCENTRATION  ALGORITHMS 


Evolution  of  passive  satellite  sensors  has  been  accompanied 
by  development  of  data  processing  techniques  and  procedures  that 
extract  information  from  radiometric  signals.  Raw  radiometer 
data  are  reduced  using  a  series  of  algorithms  that  perform 
calibration  and  analytic  functions.  Noise  is  removed,  radiome- 
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trie  corrections  are  made,  and  brightness  temperatures  are 
computed.  The  products  of  these  algorithms  typically  are  values 
of  ice  concentration,  commonly  total  ice  concentration  but 
sometimes  multi-year  and  firstyear  ice  concentrations  as  well. 
Generally  the  information  is  presented  in  map  form.  Algorithms 
exist  to  treat  both  single-channel  ESMR  data  (Gloersen,  et  al . , 
1974,  1978;  Zwally,  et  al  .  ,  1979;  Carsey,  1982  ;  Comiso  and 
Zwally,  1982)  and  multi-channel  SMMR  data  (Gloersen,  et  al . , 
1981a,  1981b;  Martin,  e t  a 1 . ,  1982;  Cavalieri,  et  al . ,  1983; 
Svendsen,  et  al . ,  1983). 

The  ability  to  estimate  ice  characteristics  from  microwave 
data  is  based  on  the  fact  that  radiation  emitted  by  ice  and  water 
is  described  by  the  Ray 1 e igh- Jeans  approximation  for  the 
intensity  of  thermal  radiation  emitted  by  a  blackbody  (Wilheit, 
1972) : 

Tb  =  ETo  (1) 

where  T_  is  the  measured  brightness  temperature  in  Kelvins,  E  is 

.  iJ 

the  emissivity  of  the  imaged  material,  and  Tq  is  the  absolute 
physical  temperature  of  the  material.  The  approximation  assumes 
that  the  radiating  medium  is  non-scattering  and  isotropic. 

Ice-laden  seas  present  complex  surfaces  in  which  several 
different  media  are  contained  within  a  single  satellite-scale 
footprint.  A  given  scene  typically  consists  of  open  water  in  the 
form  of  leads,  polynyas,  or  open  seas,  first-year  ice  and  new  ice 
of  various  thicknesses,  and  floes  composed  of  thick,  low-salinity 
multi-year  ice.  Each  of  these  materials  possesses  different 
physical  characteristics.  The  assumption  of  isotropy  thus  is 
violated.  Atmospheric  moisture,  atmospheric  emission,  and 
emission  from  free  space  further  contribute  to  total  brightness 
temperatures  sensed  by  a  satellite-borne  imager  and  compound  the 
problem  (Fig.  1).  Equation  1  must  therefore  be  rewritten  to 
account  for  each  medium  imaged  in  a  given  footprints 


tb  - 


E  T  C 
m  om  m 


+  E^T 

f  of  f 


+  Tn  C 
Bw  w 


+  T 


Ba 


( multi-year  component ) 

( first-year  component ) 
(open  water  component) 

( atmospheric  contribution ) 


(2a) 


where  C  is  the  fraction  of  the  footprint  that  a  particular 
component  encompasses  such  that: 
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(2b) 


The  total  brightness  temperature  sensed  thus  is  a  linear 
combination  of  the  relative  proportions  of  each  component  in  the 
footprint . 

Atmospheric  contributions  are  functions  of  absorption  of 
radiation  emitted  from  the  surface,  largely  by  water  vapor  (at¬ 
mospheric  opacity) ,  and  emission  of  radiation  from  the  atmosphere 
itself.  Atmospheric  effects  typically  are  small  at  the  poles  due 
to  low  humidity  conditions  characteristic  of  the  region  (Wilheit, 
et  al.,  1972).  This  is  especially  true  during  winter  months.  At 
lower  latitudes  in  the  vicinity  of  marginal  ice  zones,  open  water 
is  extensive  and  warmer  physical  temperatures  increase  the 
capacity  of  the  atmosphere  to  hold  water.  Here,  atmospheric 
attenuation  is  greater  and  its  contribution  to  brightness  temper¬ 
ature  limits  the  sensitivity  of  uncorrected  brightness  tempera¬ 
ture  measurements  to  a  larger  degree.  Variation  in  atmospheric 
water  with  latitude,  season,  weather,  and  with  the  location  of 
the  marginal  ice  zone  means  that  the  effects  of  atmosphere  can 
not  be  described  adequately  by  a  single  correction  factor  applied 
consistently,  or  even  regionally,  to  brightness  temperature 
measurements . 


FIGURE  1.  Sources  of  microwave  radiation  received  by  satellite 
sensors  (Svendsen,  et  al . ,  1983).  Brightness  temperature  is 

a  function  of  four  components:  (1)  emitted  brightness  from 
the  surface,  (2)  radiation  from  space  reflected  by  the 
surface,  (3)  upwelling  radiation  from  the  atmosphere,  and 
(4)  downwelling  atmospheric  radiation  reflected  from  the 
surface . 
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Computation  of  the  atmospheric  correction  from  radiometric 
data  requires  that  multi-channel  instrumentsObe  used.  Insuffi¬ 
cient  degrees  of  freedom  are  present  in  single-channel  data  to 
resolve  the  atmospheric  component.  Atmospheric  effects  commonly 
are  assumed  to  be  negligible  in  ESMR  ice  algorithms;  the  loss  of 
precision  that  results  is  included  in  the  overall  error  of  the 
sensor  and  the  algorithm.  With  multi-channel  data,  models  of  at¬ 
mospheric  effects  independent  of  surface  effects  can  be  developed 
(Svendsen,  et  al .  ,  1983).  The  potential  accuracy  of  algorithms 
that  utilize  multi-channel  data  thus  is  enhanced. 


ESMR  Algorithms 


The  single-channel  character  of  ESMR  data  places  limitations 
on  the  amount  of  information  it  contains.  Real-time  algorithms 
that  utilize  ESMR  data  generally  are  unable  to  discriminate 
between  first-year  and  multi-year  ice.  Delayed  analysis  of 
seasonal  data  does,  however,  yield  multi-year  ice  concentration 
measurements  (Carsey,  1982).  The  general  operational  utility  of 
ESMR  sensors  lies  in  their  proven  ability  to  map  the  ice  edge  and 
to  estimate  first-year  ice  concentrations  in  the  absence  of 
multi-year  ice.  The  simplicity  of  non-iterative,  real-time 
models  permits  reasonably  accurate  ice  forecasts  to  be  made  with 
only  modest  computational  capability. 


Real-time  Algorithms 

If  multi-year  ice  is  assumed  to  be  absent  from  the  imaged 
scene  (a  valid  assumption  in  Antarctic  waters  and  in  the  Bering 
Sea)  and  if  the  atmospheric  contribution  to  brightness  tempera¬ 
tures  sensed  by  the  satellite  is  assumed  to  be  negligible,  then 
Equation  2  can  be  simplified  to  include  only  the  first-year  and 
open  water  components : 


tb  - 


E,T  _C  _ 
f  of  f 


T  n  C 
Bw  w 


(3a) 


With  multi-year  ice  absent 
is  assumed  to  be  a  linear 
first-year  ice,  and  C^,the 


the  brightness  temperature 
combination  of  C^,  the  area 
area  covered  by  open  water: 


of  a  scene 
covered  by 


C,  +  C 
f  w 


=  1 


(3b) 


The  concentration  of  first-year  ice 
system,  then,  is  given  by: 


in  such  a  two  component 
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135 


(3c) 


Cf  = 


TB  " 


E  j-T  -  135 

f  os 


where  T  is  the  measured  brightness  temperature  of  the  scene  as 
sensed  by  the  satellite  radiometer,  135  is  the  brightness  temper¬ 
ature  of  open  water  in  Kelvins,  and  TQg  is  the  physical  tempera¬ 
ture  of  the  ice  surface  (Zwally,  et  al . ,  1979;  Comiso  and  Zwally, 
1982).  Conversely,  the  proportion  of  the  footprint  that  consists 
of  open  water  is  given  by  Gloersen,  et  al .  (1974)  as: 


C 

w 


T  D  -  Ej-T 
B  f  os 


135  -  E..T 

f  os 


(3d) 


Equations  3c  and  3d  form  the  basis  of  non-iterative  algorithms 
devised  to  produce  maps  of  first-year  ice  concentration  from  ESMR 
brightness  temperatures. 

The  accuracy  of  values  used  to  initialize  the  algorithm, 
specifically  the  emissivity  of  first-year  ice  (E^)  and  physical 
surface  temperature  (T  ),  also  limit  the  accuracy  of  predictions 
based  on  the  algorithm.  Furthermore,  where  more  than  one  ice 
type  is  present,  uncertainty  in  the  relative  proportions  of 
first-year  ice,  new  ice  types,  and  multi-year  ice  introduce 
errors.  Early  attempts  to  deduce  ice  concentration  in  the  Arctic 
basin  from  single-channel  ESMR  data  typically  failed  to  account 
for  new  ice.  Actual  concentrations  of  ice  in  the  central  basin 
typically  were  higher  than  those  predicted  from  Equation  3c  as  a 
result  of  the  influence  of  radiometr i cal ly  colder  new  ice  and 
multi-year  ice.  Because  ESMR' s  accuracy  is  limited  by  the 
presence  of  multiple  ice  types  within  the  satelite  footprint,  its 
utility  is  highest  in  the  Antacrtic  where  the  multi-year  concen¬ 
tration  within  the  pack  typically  is  small  and  the  location  of 
multi-year  ice  and  shelf  ice  can  be  predicted  accurately 
(Gloersen,  et  al . ,  1981b) . 

Overall  accuracy  for  the  model  is  estimated  to  be  +6%  for 
ice  concentration  greater  than  35%  (Gloersen,  et  al .  ,  1974)  and 
+10%  or  better  overall  (Gloersen  and  Salomonson,  1975).  These 
estimates  are  not  documented  by  ground  truth  data.  Carsey  (1982) 
sets  the  lower  limit  of  ice  concentration  that  can  be  detected 
accurately  in  the  marginal  ice  zone  from  ESMR  data  at  15%  when 
weather  conditions  are  unknown.  Ocean  roughness,  foam,  high  con¬ 
centrations  of  atmospheric  water  vapor,  and  suspended  water 
droplets  associated  with  the  occurrence  of  severe  weather 
effectively  increase  the  brightness  temperature  of  open  water  to 
that  of  a  calm  sea  with  15%  ice  concentration.  Reliable  concen- 
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tration  estimates  below  15%  total  sea  ice  thus  can  not  be  made 
unless  accurate  weather  information  is  available. 


Delay-time  Algorithms 

Carsey  (1982)  presents  a  different  approach  to  analysing 
ESMR  data  that  yields  estimates  of  the  fixed  multi-year  ice  con¬ 
centration  within  a  frozen  pack  during  a  given  winter.  Estimates 
of  multi-year  ice  obtained  from  the  algorithm  are  determined  by 
measuring  the  total  amount  of  open  water  present  at  the  end  of 
summer  melt.  Ice  that  survives  summer  melt  becomes  multi-year 
ice  upon  freeze-up.  By  determining  the  maximum  percentage  of 
open  water  at  late  summer  within  a  given  area  the  percentage  of 
winter  multi-year  ice  present  within  that  same  area  also  is 
known.  In  order  to  determine  the  maximum  open  water  area,  analy¬ 
sis  of  ESMR  brightness  temperatures  over  the  entire  freeze  period 
is  required.  Instantaneous,  real-time  estimates  of  multi-year 
ice  concentration  can  not  be  obtained  from  ESMR  data  for  this 
reason . 

Carsey  divided  the  Arctic  Basin  into  quadrants  200  km  on  a 
side.  ESMR  brightness  temperature  values  for  pixels  within  each 
quadrant  were  averaged  on  each  day  of  the  90  day  period  between 
August  12  and  November  10,  1974.  During  this  period  both  maximum 
open  water  at  the  end  of  summer  and  complete  freeze-up  at  the 
close  of  the  Arctic  autumn  occur.  Daily  brightness  temperatures 
then  were  plotted  as  a  function  of  time  for  the  entire  period  for 
each  quadrant.  The  point  at  which  maximum  open  water  occurs  is 
correlative  with  the  minimum  brightness  temperature  observed  over 
this  interval  for  each  grid  element.  The  point  at  which  all  open 
water  is  frozen  (total  freeze-up)  is  marked  by  the  maximum 
brightness  temperature  observed.  The  general  pattern  displayed 
by  each  grid  element  consists  of  a  minimum  brightness  temperature 
(maximum  open  water)  during  the  period  August  27  to  September  17. 
This  is  followed  by  a  gradual  increase  to  a  maximum  T  at  the  end 

D 

of  October.  Attainment  of  this  maximum  is  interpreted  as  an 
indication  that  freezing  of  open  water  is  complete.  If  the 
proportion  of  multi-year  ice  within  each  grid  element  is  assumed 
to  be  constant  throughout  the  freeze-up  period,  differences  in 
the  minimum  and  maximum  brightness  temperatures  observed  between 
grid  quadrants  is  a  function  of  the  relative  proportion  of 
multi-year  ice  within  each  quadrant  (Carsey,  1982).  Carsey  thus 
computes  the  concentration  of  multi-year  ice  (C  )  and  the  emis- 
sivity  of  multi-year  ice  (E  )  within  a  given  grid  element  as 
functions  of  minimum  and  maximum  brightness  temperature: 
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where  Tom(t2^  ^ s  t^ie  minlmum  physical  surface  temperature  reached 
in  a  given  quadrant  (the  physical  temperature  on  the  day  of 
maximum  open  water)  and  T  m(t^)  is  the  maximum  physical  surface 
temperature  reached  during  the  summer-winter  transition  (the 
physical  temperature  on  the  first  day  of  complete  ice  cover).  D 
is  a  correction  for  atmospheric  radiance. 

Carsey's  equations  provide  delayed  estimates  of  multi-year 
ice  concentration  for  large  (200  km)  quadrants  of  the  Arctic 
basin  on  the  basis  of  single-channel  data  taken  over  the  period 
from  maximum  open  water  (middle  to  late  August)  to  freeze-up 
(middle  to  late  October).  Neither  instantaneous  estimates  of 
multi-year  concentration,  nor  estimates  of  multi-year  concentra¬ 
tions  during  late  spring  and  early  summer  can  be  obtained 
directly  from  single-channel  data  (ESMR)  using  equations  4a  and 
4b.  However,  the  amount  of  summer  multi-year  ice  present  in  a 
given  year  probably  can  be  estimated  from  the  previous  winter's 
multi-year  cover.  If  the  concentration  of  multi-year  ice 
remaining  at  the  end  of  the  summer  is  calculated,  and  if  the  rate 
of  change  from  spring  to  fall  is  assumed  to  be  linear  (or  some 
well-behaved  function  of  climate),  then  the  multi-year  concentra¬ 
tion  in  a  specific  quadrant  can  be  estimated  after  accounting  for 
drift  of  ice  from  the  previous  year's  quadrants. 

Numeric  estimates  of  the  accuracy  of  multi-year  ice  concen¬ 
trations  computed  from  the  algorithm  are  not  provided.  Instead, 
Carsey  (1982)  makes  a  qualitative  assessment  of  potential  sources 
of  error  involved  in  the  analysis.  Variance  in  the  emissivity  of 
multi-year  ice  is  assessed  as  the  factor  that  introduces  the 
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greatest  variation.  Conditions  at  the  ice  surface  change  rapidly 
during  summer  melt  and  autumn  freeze-up  of  the  pack.  Occurrence 
of  rain  and  wet  snow  in  conjunction  with  diurnal  freezing  and 
thawing  of  water  in  the  top  centimeters  of  ice  cause  dramatic 
fluctuations  in  emissivity.  Although  radical  changes  in  bright¬ 
ness  temperature  that  result  are  well  documented  (Carsey,  1982; 
Ketchum,  1982),  relationships  between  imaged  brightness  tempera¬ 
tures  and  specific  surface  conditions  that  they  indicate  are 
poorly  understood.  The  existence  of  surface  conditions  that  give 
deceptively  cold  temperatures  near  the  time  of  maximum  open  water 
or  deceptively  warm  temperatures  when  freezing  is  nearly  complete 

would  cause  values  of  T  (t^)  and  T  (t~)  to  be  chosen 

om  2  om  3 

incorrectly . 

A  second  source  of  error  is  drift  of  ice  through  each 
quadrant  over  the  three  month  period  of  freeze-up.  The  large 
grid  size  (200  km  square)  was  chosen  such  that  areal  extent  of 
ice  that  remains  within  a  given  quadrant  for  the  entire  period 
would  be  large  compared  the  amount  of  ice  that  enters  and  leaves 
the  quadrant  as  a  result  of  drift.  At  the  same  time,  the  grid 
size  could  not  be  so  large  that  several  different  ice  regimes 
would  be  encompassed  within  each  quadrant.  Patterns  in  the 
fluctuation  of  brightness  temperatures  observed  within  individual 
quadrants  are  similar  to  those  observed  in  adjacent  quadrants. 
This  suggests  that  similar  ice  conditions  prevail  over  regions  of 
areal  extent  larger  than  the  size  of  each  quadrant  (Carsey, 
1982)..  Furthermore,  changes  in  ice  character  from  one  region  to 
another  are  transitional  rather  than  abrupt,  suggesting  that  the 
algorithm  will  remain  robust  even  if  abnormally  high  drift  rates 
move  large  quantities  of  ice  through  adjacent  quadrants. 


SMMR  Algorithms 

Algorithms  that  employ  SMMR  multi-channel  data  are  based  on 
the  premise  that  emissivities  of  different  imaged  media  (open 
water,  first-year  and  multi-year  ice,  atmospheric  water),  though 
similar  at  some  microwave  frequencies,  are  different  at  other 
frequencies.  Thus,  although  multi-year  ice  can  not  be  adequately 
discerned  from  first-year  ice  at  the  ESMR  frequency  (19.35  GHz), 
first-year  and  multi-year  ice  are  characterized  by  different 
emissivities  at  other  frequencies.  The  multi-channel  capability 
of  SMMR  instruments  enhance  the  accuracy  of  ice  concentration 
predictions  significantly.  Effects  of  multi-year  ice  and  atmos¬ 
pheric  radiance  and  opacity  can  be  accounted  for  and  their 
influence  on  predicted  total  ice  concentrations  mediated  or 
corrected.  Severe  problems  still  exist,  however,  in  transitional 
seasons  when  emissivities  of  wet  and  freezing  surfaces  are  in 
constant  flux. 
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Non-iterative  Algorithms 

Martin,  et  al .  (1982)  and  Cavalieri,  et  al .  (1983)  present  a 

two-channel,  non-iterative  algorithm  that  derives  estimates  of 
first-year  ice  concentration  from  radiation  at  the  37  GHz  fre¬ 
quency.  The  method  uses  polarization  differences  in  emissivity 
between  open  water  and  first-year  ice  to  map  the  concentration  of 
first-year  ice.  Ice  concentration  is  given  by: 

1 

C£  =  -  (5a) 

(1-R) 
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is  the  computed  concentration  of  first-year  ice  in  the  satel¬ 
lite  field  of  view,  ^ , ,  x  and  ^ and  ,,  N  and 

Bf(horz)  Bf  (vert)  Bw(hor.z) 

,  ,  x  are  known  brightness  temperatures  for  consolidated 

Bw(vert)  .  ^  .  *  , 

first-year  ice  (no  open  water)  and  water  in  horizontal  and 

vertical  polarizations,  and  T_,,  x  and  T  ,  ^  x  are  the 

B(horz)  B(ver.t) 

measured  brightness  temperatures  of  the  satellite  field  of  view 
for  the  horizontal  and  vertical  channels.  The  model  assumes  that 
atmospheric  contributions  to  imaged  brightness  temperatures  are 
negligible . 

Martin,  et  al .  (1982)  estimate  possible  errors  of  up  to  24% 

for  ice  concentration  predictions  using  this  model.  The  24% 
figure  is  based  on  ground  truth  data  acquired  during  field  exper¬ 
iments  in  the  Bering  Sea.  The  algorithm  assumes  that  only 
first-year  ice  is  present  (a  valid  assumption  for  the  Bering 
Sea).  However,  the  algorithm  does  not  discriminate  thin 
first-year  ice  from  open  water.  This  is  the  largest  source  of 
error.  Lower  concentrations  than  actually  are  present  are  given 
for  areas  with  significant  concentrations  of  thin  first-year  ice. 
Furthermore,  the  overall  accuracy  of  the  algorithm  is  dependent 
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on  the  emissivity  of  imaged  ice  types  being  close  to  the  tuned 
values  of  T  ^  and  T„w  used  in  the  algorithm.  Factors  such  as 
salinity,  surface  roughness,  and  surface  temperature  are  judged 
to  be  less  important  than  the  presence  of  ice  of  a  single  radio- 
metric  type  (Martin,  et  al . ,  1982).  Results  obtained  compare 
favorably  with  ice  edge  and  concentration  maps  produced 
independently  from  Tiros  data  and  field  data,  although  estimates 
near  landfall  are  not  as  accurate. 

NORSEX  Algorithm 

Svendsen,  et  al .  (1983)  present  an  iterative  algorithm  used 
during  NORSEX  that  predicts  both  multi-year  and  first-year  ice 
concentration  from  two  channels  of  vertically  polarized  data. 
The  37  GHz  SMMR  channel  is  used  in  conjunction  with  either  the  18 
Ghz  or  the  10.7  GHz  channel,  depending  on  the  accuracy  and 
spatial  resolution  desired.  Values  needed  for  atmospheric  cor¬ 
rections  are  extracted  from  appropriate  SMMR  channels  and  applied 
to  the  total  brightness  temperature  sensed  by  the  satellite 
before  contributions  of  the  ice-water  surface  are  determined. 

Raw  brightness  temperatures  measured  by  the  satellite  first 
are  calibrated  with  respect  to  cold  and  warm  radiators  over  areas 
of  open  ocean  and  first-year  ice  for  which  brightness  tempera¬ 
tures  are  known.  Prelaunch  calibration  data  and  polarization 
demixing  procedures  are  used  as  well.  Accuracy  of  the  NORSEX  al¬ 
gorithm  depends  in  large  part  on  the  accuracy  of  these  calibrated 
values.  Next,  an  atmospheric  correction  based  on  work  by 
Gloersen  and  Barath  (1977,  eq.  1)  is  made: 


Bs 


T0  -  2dT  t  +  dT  t  -  T 
B  a  a  a  a_ s  p 


i  -  t  -  bait  -  t  2)  -  b  -2£ 

a  a  a 

T 

a 


(6a) 


where  Tfis  is  the  corrected  brightness  temperature  contributed  by 
the  imaged  ice-water  surface  (atmospheric  effects  removed),  T  is 
the  brightness  temperature  measured  by  the  satellite  (surface 
brightness  temperature  plus  atmospheric  effects),  t  is  the  total 
atmospheric  opacity  (optical  depth  from  Reeves,  191%),  Tg  is  the 
temperature  of  free  space,  approximately  2.7  K,  and  b  Ifid  d  are 
constants  given  values  of  0.95  and  0.90  respectively.  The  term 
dT^  is  the  weighted  average  atmospheric  temperature  in  the  lower 
troposphere . 
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With  an  approximation  of  atmospheric  effects  removed,  equa¬ 
tions  2a  and  2b  are  solved  for  a  three  component  system  of 
multi-year  ice,  first-year  ice,  and  open  water  using  two  channels 
of  data: 


T0  =  C  E  T  . 
B  m  m  01 


+  C^ET  . 
f  f  01 


+  C  E  272 
w  w 


(6b) 


To.,  the  effective  physical  temperature  of  the  imaged  ice 
surface,  is  a  function  of  air  and  ocean  temperature  and  is 
approximated  by: 

T  .  =  aT  +  (1-a) 272  (6c) 

oi  oa 

where  T  is  the  physical  air  temperature  at  the  surface,  272  is 

the  physical  ocean  temperature,  and  a  is  estimated  from  field 

data  to  be  approximately  0.40.  Initial  values  for  atmosphere 

temperature  (Equation  6c,  Tq^)  are  derived  from  the  best 

available  surface  source  (ice  camps,  airborne  PRT-5  and  PRT-6 

instruments,  and  automated  ice  buoys  drifting  within  the  pack). 

Brightness  temperatures  measured  by  two  SMMR  channels  are 

used  in  Equation  6b  to  yield  two  equations  with  two  unknowns. 

These  are  solved  to  give  initial  concentrations  of  multi-year  and 

first-year  ice  (C^  and  C^).  Next,  refined  values  for  optical 

opacity  (t  )  and  atmospheric  brightness  temperature  (Tga)  are 

determined.  is  found  from  the  initial  values  of  and  C  by 

Ba  f  m  J 

interpolating  between  a  brightness  temperature  for  100%  open 

water  and  the  measured  brightness  temperature  for  100%  ice  con¬ 
centration.  The  refined  t  is  found  using  this  new  value  of  T 

t  .  oa 

b  are  solved  using  these  new  estimates 

of  T  and  t  to  produce  final  values  of  C-  and  C  . 

Ba  a  ,  ,1  m , 

Based  on  field  data  acquired  during  NORSEX  in  a  test  area 

between  Greenland  and  Svalbard,  the  NORSEX  algorithm  estimates 
total  ice  concentration  to  +3%  and  the  multi-year  concentration 
to  +^L0%  under  ideal  conditions  (Svendsen,  e t  a  1  .  ,  1983  ). 

Comparison  of  measurements  of  total  ice  concentration  derived 
from  aircraft  ground  truth  data  with  concentrations  derived  with 
the  NORSEX  algorithm  show  a  mean  difference  of  0.5^+5%  using  the 
10/3  7  GHz  SMMR  combination.  Absolute  accuracy  of  total  concen¬ 
tration  derived  from  aircraft  is  +3%.  Accuracy  of  -3.5+2%  was 
achieved  using  the  18/37  GHz  combination,  meaning  that  the  SMMR 
algorithm  estimates  concentration  an  average  of  3.5%  too  low. 
Accuracy  of  multi-year  ice  concentration  estimates  are  -4.0+6% 
for  the  10/37  GHz  combination,  and  -8.0+6%  for  the  18/37  GHz 
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combination.  Absolute  aircraft  accuracy  for  predicting 
multi-year  concentration  is  jyLO%.  The  presence  of  very  thin 
first-year  ice  or  meltponds  on  multi-year  ice  force  larger 
errors.  Wet  snow  on  multi-year  ice,  a  common  phenomenon  during 
summer  and  autumn  seasons,  causes  multi-year  ice  to  appear 
warmer,  like  first-year  ice.  Spatial  resolution  of  90  km  is 
achieved  when  the  10  GHz  and  37  GHz  channels  are  combined;  60  km 
resolution  is  achieved  using  the  18  GHz  and  37  GHz  channels. 

Gloersen  Algorithm 

Gloersen,  et  al .  (1981a,  1981b)  present  an  algorithm  for 
determining  values  of  four  variables  from  eight  channels  of  SMMR 
satellite  data.  Total  ice  concentration  (C),  multi-year  concen¬ 
tration  (F),  sea  ice  temperature  (T  ^),  and  atmospheric  liquid 
water  (L)  are  derived  using  the  vertically  and  horizontally 
polarized  brightness  temperatures  of  the  4.6,  2.8,  1.4,  and  0.8 
cm  SMMR  channels.  Brightness  temperature  is  assumed  to  be  the 
sum  of  four  discrete  sources: 


Brightness 

Temperature 


(Surface  Emissivity) 
x  (Phys.  Surface  Temp)  + 
x  (Atmos.  Absorption) 


(Upwelling 
Atmos .  Rad . ) 


(Refl.  Downwelling 
Atmos.  Radiation)  x 
(Atmos.  Absorption) 


(Refl.  Space  Rad.)  x 
(2  x  Atmos.  Absorp. ) 


Or,  in  algebraic  terms: 


Atmospheric 
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Terms 
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T  is  the  brightness  temperature  sensed  by  the  satellite  through 

channel  n,  E^n  is  the  emissivity  of  the  ocean  at  the  frequency  of 

channel  n,  E _  and  E  are  the  emissivities  of  first-year  and 
.  .  f  .  m 

multi-year  ice  at  this  same  frequency,  T  .  is  the  physical  temp¬ 
erature  of  the  ice  surface,  Tqs  is  the  physical  temperature  of 
open  space,  X  is  the  thin  first-year  ice  fraction  of  C,  and 

are  physical  atmospheric  temperatures  and  t  is  the  atmospheric 
opacity . 

Gloersen,  et  al .  (1981a)  assume  that  for  polar  regions 

and  T ^  are  approximately  equal  to  a  weighted  atmospheric  tempera¬ 
ture  of  the  lower  troposphere  which  is  approximately  equal  to  the 
temperature  of  the  ice  surface.  On  the  basis  of  this  assumption 
and  others  that  concern  the  relatively  small  opacity  of  the  polar 
atmosphere.  Equation  7a  is  converted  to: 


Bn 
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wn  f  oi 
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where  271  is  the  physical  ocean  temperature  in  Kelvins,  E^  is 
the  difference  between  the  emissivity  of  first-year  ice  and 
first-year  thin  ice,  and  1  is  the  wavelength  of  channel  n. 

Initially,  Gloersen,  et  al .  (1981a)  sought  to  solve  Equation 

7c  simultaneously  for  six  unknowns:  the  four  given  above  (C,  F, 
T  .  ,  and  L)  plus  the  thin  first-year  ice  fraction  of  C  (X)  and 

OI 

the  difference  between  thick  and  thin  ice  temperature.  However, 

high  correlations  between  brightness  temperatures  recorded  for 

adjacent  channels  produced  unacceptably  high  noise  factors  for 

many  parameters.  The  revised  algorithm  solves  iteratively  for  C, 

T  .  ,  L,  and  F,  T  L,  and  F  are  solved  at  a  resolution  of  156 
oi  oi 

km;  C  is  solved  at  a  scale  of  25  km.  Pre-launch  values  for  E  , 

o 

#  Em  #  ant  Eft  are  use<^  to  initialize  Equation  7c,  The  algo¬ 
rithm  shown  in  Figure  2  is  followed  to  compute  the  above 
variables.  Convergence  typically  occurs  in  approximately  five 
iterations  (Gloersen,  et  al , ,  1981a), 

Gloersen,  e t  a  1  ,  (1981a)  claim  to  achieve  some  rather 

remarkable  accuracies  with  their  algorithm.  Total  ice  concentra¬ 
tion  at  a  156  km  scale  is  estimated  to  be  +4%  after  the  first 
iteration,  and  as  good  as  +1%  after  several  iterations.  Total 
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FIGURE  2.  Algorithm  for  determining  total  ice  concentration 
(C),  sea  ice  temperature  (T  .),  atmospheric  liquid  water 
(L),  and  multi-year  ice  concentration  (F)  from  SMMR  data 
(Gloersen,  et  al . ,  1981a) .  Brightness  temperatures  txat  are 
input  and  the  size  of  their  footprints  are  given  in  boxes  to 
the  left.  Parameters  shown  to  the  right  are  derived  at  re¬ 
solutions  indicated  from  these  inputs  using  Equations  7a 
through  7c  (see  text). 
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concentrations  at  the  100,  50,  and  25  km  scales  after  iterating 
are  estimated  to  be  +2,  3,  and  5%.  Physical  ice  temperature  is 
estimated  to  an  accuracy  of  9  K  initially,  and  3  K  ultimately  at 
the  156  km  scale.  Atmospheric  liquid  water  is  estimated  to 
443.015  cm  (+0.011  cm  with  iteration)  and  the  multi-year  fraction 
of  total  ice  concentration  to  +24%  ( jd.9%  after  iterating)  ,  also 
at  the  156  km  scale.  Ground  truth  data  to  support  these  accuracy 
estimates  are  not  presented  by  Gloersen,  etal.  (1981a,  1981b). 


SSM/I  IMPACT  ON  ICE  ANALYSIS 

SMMR  and  ESMR  products  provide  indications  of  the  minimum 
accuracy  and  versatility  that  can  be  expected  from  SSM/l  data. 
They  also  expose  weaknesses  in  currently  used  analysis  strategies 
that  should  be  strengthed  before  the  next  generation  of  satellite 
imagery  is  available.  Review  of  past  work  is  instructive  when 
viewed  in  this  light.  If  improved  estimates  of  ice  pack  charac¬ 
teristics  are  to  be  extracted  from  SSM/l  data,  then  improved  data 
analysis  techniques  must  be  developed  and  implemented.  Previous 
applications  for  satellite-borne  microwave  imagery  fall  into  four 
general  categories:  (1)  ice  concentration,  (2)  location  of  the 
ice  edge  and  characteristics  of  the  marginal  ice  zone,  (3) 
location  of  polynyas,  large  leads,  and  open  water,  and  (4)  esti¬ 
mates  of  heat  exchange  in  polar  regions.  The  state  of  each  is 
assessed  below. 


Ice  Concentration 

Estimates  of  ice  concentration  based  on  ESMR  and  SMMR  data 
are  derived  from  algorithms  described  above.  SMMR  algorithms 
should  be  applicable  to  SSM/l  data  with  alterations  to  accomodate 
differences  in  frequency.  The  SMMR  sensor  suite  includes  radio¬ 
meters  at  10,  18,  21,  and  37  GHz;  SSM/l  will  carry  radiometers  at 
19,  22,  37,  and  85  GHz.  Emissivities  of  multi-year  and 
first-year  ice  at  85  GHz  are  similar  to  those  at  37  GHz,  the  high 
frequency  channel  aboard  SMMR  instruments  (Fig.  3)  (Troy,  et  al . , 
1981).  The  improved  spatial  resolution  that  85  GHz  offers  (14.8 
km  vs.  33  km)  will  enhance  the  quality  of  SSM/l  products  that  are 
based  primarily  on  this  sensor:  location  of  the  ice  edge  and 
determination  of  total  ice  concentration  in  the  absence  of 
multi-year  ice.  Therefore,  data  from  the  85.5  GHz  channel  should 
be  utilized  in  computing  these  variables.  Discrimination  between 
thick  first-year  ice  and  new  ice  or  thin  first-year  ice,  however, 
is  ambiguous  at  85  GHz  (Troy,  et  al  .  ,  1981)  .  Use  of  high  fre¬ 
quency  channels  in  conjunction  with  lower  resolution,  lower  fre¬ 
quency  channels  will  be  necessary  to  discern  first-year  ice 
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thicknesses  and  to  compute  multi-year  ice  concentrations. 
Spatial  resolution  for  these  estimates  could  be  as  poor  as  55  km 
if  data  from  the  19  GHz  channel  are  utilized. 

Ice  concentration  estimates  derived  from  SSM/l  data  will  be 
comparable  in  accuracy  to  those  attained  using  SMMR  imagery. 
Svendsen,  et  al .  (1983  )  present  the  best  documented  SMMR  algo¬ 
rithm  estimates  (NORSEX)  .  Total  concentration  is  estimated  to 
+3%  and  multi-year  ice  concentrations  are  estimated  to  accuracies 
of  +^10%  using  the  10  GHz  and  37  GHz  vertical  channels  combined. 
These  errors  are  based  on  ground  truth  data  and  represent  figures 
for  optimum  conditions.  Where  very  thin  ice  is  present  or  where 
meltponds  occur,  errors  will  be  larger  in  the  absence  of  ground 
truth  data.  Furthermore,  multi-year  ice  covered  by  wet  snow  will 
appear  as  first-year  ice  and  increase  the  error  of  multi-year 
concentration  estimates  (Svendsen,  et  al.,  1983). 


FREQUENCY  (GHz) 


FIGURE  3.  Emissivity  as  a  function  of  frequency  between  14  and 
90  GHz  for  four  ice  types  (from  Troy,  et  al . ,  1981). 

Vertical  bars  mark  one  standard  deviation  from  the  plotted 
mean . 
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In  general,  the  accuracy  of  concentration  estimates  will  not 
be  improved  by  the  higher-resolution  sensor  suite  aboard  SSM/l. 
Although  exclusive  use  of  the  85  GHz  channel  to  estimate  total 
ice  concentration  will  reduce  pixel  size  (area)  by  approximately 
75%,  errors  in  approximating  actual  concentrations  within  each 
cell  will  be  unaffected.  The  accuracy  of  a  given  estimate  is 
largely  a  function  of  the  algorithm  by  which  it  is  produced. 
Only  by  improving  analysis  procedures  will  the  accuracy  of  esti¬ 
mates  be  increased. 

The  ultimate  accuracy  of  any  ice  classification  algorithm  is 
determined  principally  by  three  factors:  (1)  the  information 
content  of  the  raw  data,  (2)  the  ability  of  the  algorithm  to 
remove  noise  from  the  raw  data,  and  (3)  the  ability  of  the  algo¬ 
rithm  to  discern  ice  types  on  the  basis  of  information  contained 
in  raw  data  from  which  noise  has  been  removed.  The  information 
content  of  data  gathered  by  a  particular  suite  of  sensors  is 
dependent  on  radiometer  frequencies  and  polarizations  incorpo¬ 
rated  into  design  of  the  instrument.  Characteristics  of  the  in¬ 
strument  in  effect  set  boundary  conditions  that  determine  the 
type  of  information  that  can  be  extracted  as  well  as  the  upper 
limits  of  the  quality  of  that  information.  Thus,  the  highest 
level  of  accuracy  that  can  be  obtained  in  predicting  a  given 
parameter  is  determined  by  innate  characteristics  of  the  sensor 
suite.  The  choice  of  sensors  is  invariant  after  the  design 
freeze  and  is  a  fixed  quantity  that  can  not  be  altered.  Improved 
accuracy  of  SSM/l  data  over  SMMR  data,  then,  will  only  be 
achieved  to  the  extent  that  noise  removal  techniques  can  be 
enhanced  and  discrimination  strategies  that  are  used  to  recognize 
specific  ice  types  are  improved. 

Noise,  as  used  here,  is  a  function  of  environmental  factors 
that  create  differences  between  the  brightness  temperature 
perceived  by  the  satellite  sensor  and  energy  actually  emitted  by 
the  ice-water  surface.  Atmospheric  opacity,  atmospheric  emis¬ 
sion,  emission  from  free  space,  and  inaccurate  calibration 
correction  procedures  are  the  primary  sources  of  error  associated 
with  noise.  Atmospheric  corrections  either  are  based  on  assumed 
values  of  atmospheric  opacity  and  emission  (Gloersen,  et  al .  , 
1981a)  or  are  determined  iteratively  from  multi-channel  data 
(Svendsen,  et  al .  ,  1983).  Some  workers  assume  that  atmospheric 
effects  are  minimal  and  fail  to  make  any  correction  at  all 
(Martin,  et  al . ,  1982;  Cavalieri,  et  al . ,  1983).  Such  algorithms 
are  neither  as  accurate  nor  as  reliable  as  algorithms  in  which 
atmospheric  effects  are  removed  completely. 

Calibration  corrections  are  applied  to  raw  data  received 
from  the  satellite  before  processing  by  ice  classification  algo¬ 
rithms  and  also  affect  accuracy.  Raw  data  received  from  the  sat¬ 
ellite  are  in  the  form  of  voltages  that  are  proportional  to 
differences  between  the  sensed  radiometric  signal  and  internal 
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reference  sources  aboard  the  satellite  (Swanson  and  Riley,  1980). 
Brightness  temperatures  are  determined  from  these  voltage  data 
using  calibration  algorithms.  Pre-launch  data  that  describe  the 
behavior  of  the  instrument  are  used  in  conjunction  with 
polarization  demixing  procedures  to  convert  received  voltages  to 
brightness  temperatures  (Gloersen  and  Barath,  1977;  Njoku,  et 
al .  ,  1980b;  Swanson  and  Riley,  1980).  Additional  corrections 
based  on  measurements  over  warm  and  cold  radiators  of  known  value 
are  used  by  some  workers  to  increase  further  the  accuracy  of 
computed  brightness  temperatures  (Svendsen,  e t  a  1 .  ,  1983  ). 
Clearly,  the  accuracy  of  these  corrections  and  calibrations  is  a 
critical  factor  in  determining  the  ultimate  accuracy  with  which 
ice  parameters  are  estimated. 

The  ability  of  the  algorithm  to  discern  ice  types  from  these 
corrected  signals  is  a  function  of  three  factors:  the  overall 
discrimination  strategy  from  which  the  algorithm  is  defined,  the 
accuracy  of  values  used  to  initialize  iterative  algorithms,  and 
fundamental  assumptions  made  regarding  the  type  of  materials 
present  in  the  imaged  field  of  view.  Existing  algorithms  fail  to 
describe  summer  and  autumn  conditions  adequately.  Fluctuations 
in  the  emissivity  of  summer  snow-covered  ice  caused  by  changing 
weather  conditions  and  diurnal  freeze-thaw  cycles  are  poorly 
understood.  Accordingly,  discrimination  strategies  are  poorly 
formulated,  assumptions  concerning  the  radiometric  appearance  of 
imaged  ice  are  violated,  and  seed  values  used  to  initialize 
certain  algorithms  are  poorly  chosen.  Significant  improvements 
in  the  overall  performance  of  SSM/l  algorithms  can  be  made  by 
formulating  strategies  designed  to  improve  predictions  of  ice 
concentration  during  summer  melt  and  autumn  freeze-up. 


Ice  Edge  and  Marginal  Ice  Zone 


A  primary  use  of  satellite  remote  sensing  in  arctic  regions 
is  location  of  the  ice  edge.  Daily  updates  derived  from 
microwave  sensors  provide  information  needed  to  keep  ice  charts 
current.  Microwave  imagery  provided  by  ESMR  and  SMMR  type  in¬ 
struments  represent  significant  improvements  over  visible  and 
infrared  imagery  in  that  they  return  images  of  the  surface  re¬ 
gardless  of  weather  and  light  conditions.  Accuracy  in  locating 
the  edge  is  determined  by  sensor  resolution  and  by  ice  and 
weather  conditions  at  the  ice  edge.  Differences  between  the 
emissivity  of  water  and  ice  make  the  delineation  of  compact  ice 
edges  reasonably  unambiguous.  Dispersed  edges  are  located  less 
easily.  The  minimum  ice  concentration  that  can  be  detected 
reliably  is  15%  in  the  absence  of  accurate  weather  data  regard¬ 
less  of  sensor  resolution  (Carsey,  1982).  An  effective  resolu¬ 
tion  of  30  km  per  square  pixel  is  achieved  by  most  microwave 
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imagers  to  date.  SSM/l's  85.5  GHz  channel  with  a  footprint  14.8 
km  square  should  enhance  our  ability  to  locate  the  edge  more 
accurately.  SMMR  accuracies  comparable  to  or  better  than  those 
achieved  with  Tiros  data  are  reported  by  several  workers  (Martin, 
et  al  .  ,  1982;  Cavalieri,  e  t  a  1  .  ,  1983  ).  Work  performed  by 

(1983)  in  conjunction  with  NORSEX  located  the 
+  1 0  km  from  SMMR  data.  Improved  accuracy  of  ap¬ 
is  expected  of  SSM/l  based  on  the  improved  resol- 
5  GHz  channel. 

of  specific  features  within  the  marginal  ice  zone 
A  change  in  brightness  temperature  of 


Svendsen,  et  al . 
edge  to  within 
proximately  50% 
ution  of  the  85. 

Resolution 

is  achieved  less  easily, 
more  than  60  K  over  a  100  km  wide  zone  is  required  to  insure  that 
geographical  averaging  across  the  30  km  square  ESMR  field  of  view 
does  not  obscure  features  within  the  marginal  pack  (Carsey, 

1982) .  As  a  consequence,  detailed  pictures  of  ice  concentration 
have  not  been  developed,  even  with  SMMR  data  (Svendsen,  et  al .  , 

1983) .  Nonetheless,  some  aspects  of  the  character  of  ice  within 

the  zone  can  be  deduced.  SMMR  sensors  show  differences  in 
response  between  compact  and  dispersed  edges  making  qualitative 
descriptions  of  ice  conditions  possible.  Some  indication  of  the 
type  of  ice  present  also  is  suggested  (Svendsen,  et  al .  ,  1983). 

Improved  resolution  offered  by  the  85.5  GHz  SSM/l  channel  should 
permit  improved  concentration  estimates  to  be  made  for  the  margi¬ 
nal  ice  zone.  However,  determination  of  ice  age  information  is 
likely  to  remain  out  of  reach  as  a  result  of  the  low  resolution 
of  lower  frequency  sensors  needed  to  discern  multi-year  ice  con¬ 
centrations  accurately . 
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Polynyas 

Several  groups  of  workers  have  studied  large  polynyas  using 
ESMR  imagery.  Zwally  and  Gloersen  (1977),  Rayner  and  Howarth 
(1979),  and  Carsey  (1980)  report  on  the  large  (2  x  10^  km^) 
recurring  polynya  in  the  Weddell  Sea  offshore  from  Antarctica. 
Carsey  (1980)  monitored  development  and  decay  of  the  polynya 
through  four  winters.  Changes  in  shape,  location,  and  size  were 
developed  from  ESMR  images.  Crawford  and  Parkinson  (1981) 
monitored  winter  changes  in  the  shape  and  size  of  North  Water 
Polynya  for  three  winter  periods.  Crane,  et  al .  (1982)  and  Crane 
(1983)  report  relationships  between  late  summer  polynyas  in  the 
Arctic  Basin  and  regional  weather  patterns.  Cavalieri,  et  al . 
(1983)  describe  the  development  of  large  lee-shore  polynas  in  the 
Bering  Sea  caused  by  the  movement  of  pack  ice  past  St.  Matthew 
and  St.  Lawrence  Islands.  This  last  study  used  SMMR  imagery  and 
demonstrated  that  the  location  of  such  areas  of  low  ice-concen¬ 
tration  can  be  documented  more  accurately  with  SMMR  imagery  than 
with  conventional  Tiros  data.  The  overall  accuracy  of  locating 
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polynya  boundaries  is  similar  to  accuracies  involved  in  locating 
ice  edges  (+10  km,  Svendsen,  et  al . ,  1983). 

Each  of  these  studies  reports  on  large  polynyas.  Smaller 
features  are  difficult  to  detect  using  ESMR  and  SMMR  sensors  due 
to  the  poor  instrument  resolution.  Many  tactically  significant 
leads  and  polynyas  are  likely  to  be  smaller  than  the  15  km 
footprint  of  the  highest  frequency  SSM/l  sensor.  Although 
averaging  across  footprints  that  contain  large  percentages  of 
open  water  will  produce  differences  between  brightness  tempera¬ 
tures  seen  in  surrounding  pixels,  noise  in  the  imaging  system 
could  easily  mask  subtle  differences  needed  to  detect  tactically 
significant  areas  of  open  water  reliably. 

Integrated  brightness  temperature  measurements  carry 
insufficient  information  to  determine  the  configuration  of  open 
water  bodies  contained  within  single  pixels  if  conventional 
ESMR-SMMR  analytic  procedures  described  above  are  used. 
Brightness  temperatures  that  indicate  50%  open  water  could 
reflect  a  single  circular  polynya  at  the  pixel  center  as  well  as 
an  anastomosing  network  of  narrow  leads,  the  total  area  of  which 
is  50%  of  the  field  of  view.  Image  enhancement  techniques  that 
utilize  multiple  images  of  the  same  feature  taken  on  successive 
orbits  conceivably  could  be  developed  to  reduce  the  clouding 
effects  of  signal  noise  and  improve  the  reliability  of 
single-pixel  information. 


Arctic  Heat  Exchange  and  Energy  Balance 

Carsey  (1982)  and  Svendsen,  et  al .  (1983)  demonstrate  the 
utility  of  using  satellite  derived  microwave  imagery  to  estimate 
heat  exchange  between  ocean  and  atmosphere  in  polar  regions. 
Changes  in  ice  concentration  are  known  to  affect  synoptic  weather 
patterns  (Ackley  and  Keliher,  1976;  Carleton,  1981;  Cavalieri  and 
Parkinson,  1981;  Crane,  et  al . ,  1982;  Crane,  1983).  Energy  flux 
drives  regional  weather  patterns  and  is  a  function  of  both  the 
areal  extent  of  insulating  ice  that  covers  warm  ocean  waters  and 
the  thickness  of  that  cover.  The  concentration  of  multi-year  ice 
present  at  a  given  time  is  used  to  determine  these  parameters. 
Both  ESMR  (Carsey,  1982)  and  SMMR  (Svendsen,  et  al.  ,  1983)  data 
have  been  used. 

Accurate  data  regarding  multi-year  ice  covers,  though 
difficult  to  obtain,  are  crucial  to  the  accuracy  of  global  models 
that  forecast  long  term  weather  patterns  (Carsey,  1982).  When 
the  total  ice  concentration  is  high,  as  it  is  during  winter 
months,  small  errors  in  the  estimate  of  total  ice  concentration 
result  in  large  errors  in  the  heat  budget  estimate.  For  example, 
a  3%  error  for  total  ice  concentration  of  85%  produces  a  heat 
budget  error  of  15%  (Maykut,  1978;  Svendsen,  et  al .  ,  1983).  If 
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the  actual  total  ice  concentration  is  95%,  an  error  in  the  SMMR 
estimate  of  3%  produces  an  error  of  50%  in  the  heat  exchange 
estimate.  Improvement  in  accuracy  of  multi-year  ice  concentra¬ 
tions  derived  from  SSM/l  algorithms  thus  would  represent  a 
significant  improvement  in  estimates  of  global  heat  balance. 


CONCLUSIONS  AND  RECOMMENDATIONS 


Imagery  provided  by  the  SSM/l  sensor  will  represent  marginal 
improvement  over  SMMR  data,  and  significant  improvement  over  ESMR 
data.  Although  resolution  of  the  high  frequency  85.5  GHz  SSM/l 
channel  is  substantially  greater  than  the  30  km  resolution 
offered  by  SMMR's  highest  frequency,  data  received  from  this 
channel  will  affect  the  accuracy  of  only  some  ice  analysis 
products.  Spatial  resolution  of  the  position  of  ice  edges  and 
the  location  of  open  water  within  the  pack  will  be  improved. 
Further,  the  cell  size  for  which  total  ice  concentration  can  be 
estimated  will  be  reduced.  The  85.5  GHz  channel  should  be  used 
to  compute  these  parameters.  However,  improvements  in  the  accu¬ 
racy  with  which  multi-year  ice  concentration,  the  concentration 
of  thin  ice  types,  and  the  character  of  ice  in  the  marginal  ice 
zone  are  predicted  must  come  from  improvements  in  ice  algorithms. 
Improved  ground  truth  information  from  which  satellite  data  can 
be  interpreted  and  verified  is  critical  if  predictive  capabili¬ 
ties  are  to  be  improved. 

In  this  light,  research  that  addresses  the  following  points 
should  be  supported  if  the  utility  of  ice  analysis  products  is  to 
be  enhanced: 

1.  Information  from  the  85.5  GHz  channel  must  be  incorpo¬ 
rated  into  existing  algorithms  to  capitalize  on  the  improved  re¬ 
solution  that  this  sensor  offers.  The  Hughes  algorithm  presently 
implemented  at  Fleet  Numerical  Oceanography  Command  (FNOC) 
facilites  does  not  utilize  data  from  this  channel. 

2.  Current  algorithms  designed  to  treat  SMMR  data  probably 
can  be  modified  to  accept  SSM/l  data,  given  the  similarity  of 
frequencies  involved.  It  might  be  equally  efficient  to  write  new 
algorithms  that  overcome  weaknesses  of  the  SMMR  algorithm  suite 
discussed  here.  The  approach  generally  taken  in  the  past  was  to 
design  one  algorithm  that  provides  estimates  for  all  parameters. 
The  possibility  of  developing  specialty  algorithms  that  yield 
better  estimates  of  specific  parameters  from  SSM/l  data  in 
conjunction  with  data  from  other  satellite  sensors  should  be 
explored . 

3.  Radiometric  signatures  of  wet  or  snow-covered  ice 
surfaces  that  occur  during  summer  and  autumn  months  are  not  well 
documented.  Brightness  temperature  data  from  which  predictive 
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4. 

on  separate  orbital  passes 
lite-borne  microwave  imagery, 
sensor  resolutions  might  be 


algorithms  can  be  written  to  discriminate  multi-year  ice  from 
first-year  ice  during  melt  and  freeze-up  are  needed. 

Processing  techniques  that  examine  common  areas  viewed 

have  not  been  applied  to  satel- 
Features  that  occur  at  scales  near 
enhanced  if  such  techniques  are 
implemented.  Conditions  within  the  marginal  ice  zone  and  tacti¬ 
cally  significant  areas  of  open  water  and  thin  ice  of  areal 
extent  smaller  than  can  now  be  resolved  might  be  identified  from 
SSM/l  images. 

5.  Ground  truth  data  that  correlate  marginal  ice  zone 
conditions  with  imagery  from  satellite-borne  microwave  instru¬ 
ments  are  needed  to  improve  estimates  of  conditions  within  the 
marginal  ice  zone.  The  possibility  of  developing  specialty  algo¬ 
rithms  designed  specifically  to  treat  data  within  the  marginal 
ice  zone  should  be  explored. 

6.  Current  algorithms  fail  to  discriminate  adequately 
between  thin  ice  types.  Ground  truth  data  that  describe  the 
transition  from  open  water  to  thick  first-year  ice  are  needed. 


; 
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